Deterministic quantum teleportation between distant atomic objects 
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Quantum teleportation is a key ingredient of quantum networks 1 1 2 1 and a building block for quantum com- 
putation (SJHI. Teleportation between distant material objects using light as the quantum information carrier has 
been a particularly exciting goal. Here we demonstrate a new element of the quantum teleportation landscape, 
the deterministic continuous variable (cv) teleportation between distant material objects. The objects are macro- 
scopic atomic ensembles at room temperature. Entanglement required for teleportation is distributed by light 
propagating from one ensemble to the other. Quantum states encoded in a collective spin state of one ensemble 
are teleported onto another ensemble using this entanglement and homodyne measurements on light. By im- 
plementing process tomography, we demonstrate that the experimental fidelity of the quantum teleportation is 
higher than that achievable by any classical process. Furthermore, we demonstrate the benefits of deterministic 
teleportation by teleporting a dynamically changing sequence of spin states from one distant object onto another. 



Quantum teleportation of discrete and continuous 161 
variables is the transfer of a quantum mechanical state with- 
out the transmission of a physical system carrying this state. 
The first experimental teleportation protocols employed light 
as the carrier of quantum states |7, 8|. Teleportation of 
atomic states over microns distances has been realized in two 
experiments using short range interactions between trapped 
ions |9|. Interspecies teleportation from light onto atoms 
has been achieved both deterministically for continuous vari- 
ables ifTOl and probabilistically for discrete variables 111 11 . 
Recently, probabilistic teleportation between two ions |[T2ll 
and two atomic ensembles 1 13] over a macroscopic distance 
has been demonstrated. While probabilistic teleportation, in 
which entanglement is distributed by photon counting I 71IT41. 
is capable of reaching distances of many km ||21[T5|, the power 
of continuous variable (cv) teleportation is that it succeeds de- 
terministically in every attempt EKTOl, that it is capable of 
teleporting complex quantum states |16|, and that it can be 
used in universal quantum computation [4J. Here, we present 
an experimental demonstration of deterministic cv teleporta- 
tion between distant macroscopic atomic ensembles, thus ex- 
tending the powerful cv teleportation IH [TOl 116 1 onto atomic 
memory states. The protocol which succeeds in every attempt 
allows us to teleport dynamically changing quantum states of 
collective atomic spins with a bandwidth of tens of Hz. 

The quantum teleportation process begins with the creation 
of a pair of entangled objects. In our experiment these two 
objects are an atomic ensemble at site B and a photonic wave 
packet generated by interaction of this ensemble with a driv- 
ing light pulse (Eig. [T^). The wave packet travels to site A, 
the location of the atomic ensemble whose state is to be tele- 
ported. This step establishes a quantum link between the two 
locations. Following the interaction of the ensemble A and 
the wave packet, the cv Bell measurement is performed on the 
transmitted light. The result of this measurement is commu- 
nicated via a classical channel to site B, where it is fed back 
via local operations on the second entangled object, i.e. the 
ensemble B, thus completing the process of teleportation. 

Continuous variable teleportation is described in the lan- 
guage of canonical operators. For atomic ensembles, fully 



spin polarized along the x-axis, those are scaled dimen- 
sionless Cartesian components of the collective spin: x = 
Jy/V< Jx > and p = Jz/V< Jx > where Jx,y,z = 
YliJx.y.z (summed over all atoms i) is the collective angu- 
lar momentum of the ensemble. Here, we employ ^^^Cs 
atoms polarized into \F = 4, mi? = 4) electronic ground 
state. The usual link between the ladder operator b for col- 
lective atomic excitations |17 | of the state ttif = 3 (Fig.[T])) 
and canonical variables is = {x — ip)/V2. Atoms are 
placed in a bias magnetic field collinear with J^, so that 
in the lab frame the observables x oc Jy and p oc Jz ro- 
tate at Larmor frequency ft (Fig. [it) according to the atomic 
Hamiltonian Ha = — ^ {x'^ ^ P) /2. Atom-light interac- 
tion is shown in Fig. [T]) and involves two scattering pro- 
cesses Hint oc /J^alsb^ — ^<^is^ + h.c. where aj^/j^ generate 
photons in the upper/lower (u ± sideband modes of the 
driving field u. Previously, this type of interaction has been 
used to create entanglement between photonic modes iTTSl , 
described by aus,«is, and between two atomic spin ensem- 
bles |[T9ll . This interaction contains both essential ingredi- 
ents of the teleportation protocol, the creation of entangle- 
ment between the atoms and photons (first term) and a beam- 
splitter type operation between atoms and photons (second 
term). For the present experiment, the ratio of the amplitudes 
of the entangling and beam-splitter parts of the Hamiltonian 
is jj^/u = 1.38. Due to the Larmor precession, it is useful 
to introduce two-mode canonical variables for light (y and g), 
given by i/c cos{rtt) + i/s s'm{Qt) ex ausC"*^^ + ai^e'^^^ + h.c. 
where ?/c,s are the cosine and sine components, with similar 
definitions for q. Encoding quantum states of atoms in high 
frequency modes ft allows for quantum noise limited perfor- 
mance of cv variables of macroscopic atomic ensembles. Pho- 
tonic variables i/c^s are measured via polarization homodyne 
detection with nearly unity quantum efficiency and are also 
quantum noise limited at frequency Q. 

The experiment (Fig.[T^) utilizes two room temperature gas 
ensembles of Na ^ 10 — 10^^ spin polarized Cesium atoms 
in glass cells with spin protecting coating as in |[171 118J |20l 
placed at a distance of 0.5m. Crucially, the macroscopic 
spins of the two ensembles here are parallel, whereas in our 
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Figure 1: Teleportation experiment, (a) The experimental layout. 
A strong driving pulse propagates first through ensemble B creating 
the mode y' , q entangled with B and then through ensemble A whose 
state is to be teleported. Bell measurement on the modes y'\ q" is 
performed by polarization homodyning. Teleportation is completed 
by classical communication of these results to B. More comments 
in the text, (b) The level scheme and relevant transitions. Classi- 
cal drive field (thick lines) and quantum fields forming the mode ^, q 
(wavy lines) are shown, (c) The time line of the experiment. I. Prepa- 
ration of the input state, II. Entanglement and Bell measurement. III. 
Feedback, IV. Read-out of the teleported state. 



previous work (TtI [191 where the goal was generation of 
entanglement between two ensembles they were antiparal- 
lel. The bias magnetic field of B ^ 0.9G corresponds to 
Vt — 322kHz. Both ensembles are initially optically pumped 
into the \F = 4,mi? = 4) state using the combination of 
852nm and 893nm circularly polarized light, thus creating a 
coherent spin state (CSS) : Var(J^) • Var(J;,) = jJ/4 with 
^ 4A^^ and {Jy) — and (J^) = 0, corresponding 
to a vacuum state with variances Var(x) = Var(p) = 1/2. 
The spin of the ensemble A to be teleported is then displaced 
with mean values {xa] and (p^) by a weak radio-frequency 
(rf) magnetic field pulse of frequency corresponding to the 
creation of coherent superposition of electronic ground states 
vfip =3, mi? = 4 (Fig.[TJ)). This input state can also be made 
dynamically evolving in time, as described below. 

The layout and time sequence for teleportation and veri- 
fication are shown in Fig. [T^,c. A ^/-polarized, 3ms long 
light pulse drives the interaction with a power of 5.6mW, 
blue detuned from the line F = 4 ^ 5 transition by 
A = 850MHz. The forward scattered mode, x-polarized 
and described by y' is entangled with the collective spin 
B and co-propagates with the drive light towards the site A. 
The interaction with the ensemble A leads to partial map- 
ping of its state onto light, and is followed by the Bell mea- 
surements on the light modes y" ^ q" of the upper and lower 
sidebands performed via polarization homodyning with the 
driving light acting as the local oscillator. The detailed the- 
oretical model of the teleportation protocol is presented in 
the Supplementary Information (SI). Near unity teleporta- 
tion fidelity can be achieved |23|, if the driving fields for A 
and B ensembles are made time-dependent. However, even 
with top hat driving pulses a sufficiently high fidelity can 
be achieved, if an optimal temporal mode for the detected 



homodyne signal is chosen. The optimal readout mode is 
Uc/s,- oc Jq COS I sm{ftt)e~'^^y{t)dt, where T is the pulse 
duration and 7 is the decay rate of the atomic state. Measure- 
ments of such observables with specific temporal mode func- 
tions can easily be conducted by electronic post-processing of 
the photocurrent. The teleportation protocol is completed by 
sending the measurement results y^^^^ _ via a classical link to 
the site B where spin rotations in they^z plane conditioned on 
these signals are performed using phase and amplitude con- 
trolled rf magnetic field pulses at frequency ft. The gain coef- 
ficient for this classical feedback is optimized for every class 
of teleported states to maximize the average fidelity. The de- 
terministic character of the homodyne process ensures success 
of the teleportation in every attempt. 

We performed teleportation of sets of CSS with varying 

< Jy,z > , corresponding to displaced vacuum or coherent 
states in quantum optical terms. Calibration of the input states 
of ensemble A and the post-teleportation read-out of the re- 
sulting state of ensemble B are done by polarization homo- 
dyning of transmitted pulses, y^^^ _ and y^^^ _, respectively, 
with the same temporal mode as for the Bell measurement. 
Using Eq[T] from the Methods Summary, the mean values of 
the input states of ensemble A are found as < xa >=< 
y^_ > Ik, < Pa >=< yf_ > / k and their variances 
as'var(xA) = (Var(7/,^_) -'4/2 - cl/2 - c%m/2)/K\ 
Var(j9A) = (Var(y^_)-cg/2-c^/2-c^m/2)//^2 with inter- 
action coefficients Cy^q^N-,'^ and k described in the Methods 
Summary. As in our previous work with room temperature 
spin ensembles, the prepared atomic input states are found 
to be very close to ideal CSS' with Var(x^) = Var(pA) = 
(1.03 ±0.03) -1/2, which confirms the validity of the read-out 
procedure and allows for teleportation process tomography. 
After each teleportation sequence, the characteristic variables 

< x^^^ >, < p^^^ >, Var(x^^^^), Var(p^^^^) of the spin state of 
the target ensemble B are found in the same way from the 
read-out of the verification pulse y^j^ _. Based on those val- 
ues the experimental fidelity is determined using a standard 
method of calculation of state overlap |23|. Optimization of 
the teleportation protocol has been performed by varying the 
drive pulse duration T, the measured temporal mode of light, 
and the gain for the classical feedback. The optimal read-out 
mode was always found with an exponential decay rate equal 
to the spin decay 7 as expected from the model. 

The quantum character of the teleportation is verified by 
comparing the fidelity of state transfer to the classical bench- 
mark teleportation fidelity. More specifically, we perform the 
teleportation process tomography 1 10, 21] using various sets 
of coherent spin states of ensemble A as input states. For such 
Gaussian states the individual state transfer fidelity is calcu- 
lated from the first two moments |23|, i.e. the mean values 
and the variances = Var(x^^^^), cr^ = Yar{p^^^). We then 
evaluate the average transfer fidelity for sets of coherent states 
with a Gaussian distribution of displacements with mean num- 
ber n =< b^b > of spin excitations | 22, 23 1. A rigorous clas- 
sical benchmark fidelity (1 + n)/(l + 2n) for such classes of 
states has been derived in |22|. Demonstration of a fidelity 
above the classical benchmark signifies the success of quan- 
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Figure 2: Teleportation fidelity, (a) The variances of the teleported 
state Var(x^^^^,p^^^^) in projection (vacuum) noise units as a function 
of gain. (b)Teleportation fideUty as a function of the mean photon 
number of the input set of states. Blue curves/points - feedback ap- 
pUed before the verification measurement, red curves/points - feed- 
back apphed after the verification measurement. Black dashed line 
represents the classical benchmark. Error bars represent one standard 
deviation. 



The grey points represent results for individual teleportation 
runs. 

The fidelity of the teleportation can be further improved by 
using time varying drive pulses |23 1 and increasing the optical 
depth of the atomic ensembles. Cv teleportation is capable of 
teleporting highly non-classical states as shown for teleporta- 
tion of light modes 1 16 1, so it can be expected that determinis- 
tic teleportation of an atomic qubit 1 14 | can be performed by 
developing the present approach. The stroboscopic teleporta- 
tion of spin dynamics can also be extended towards a true con- 
tinuous in time teleportation paving the way to teleportation of 
quantum dynamics and simulations of the interaction between 
two distant objects which do not interact directly |24|. From 
the quantum sensing perspective teleportation of spin dynam- 
ics allows for performing measurements at remote locations, 
spatially separated from the environment where the actual dy- 
namics take place. 

Methods Summary 



tum teleportation and is equivalent to the ability of the telepor- 
tation channel to transfer entangled states. The Gaussian char- 
acter of the atomic spin states has been experimentally veri- 
fied, and experimental runs for different CSS input states with 
displacements of 0, 5, 25, 160 in canonical units and phases 
0,7r/4,7r/2inx,p space and different gains g have been per- 
formed. For every input state, 10.000-20.000 teleportations 
can be performed within a time slot of a few minutes with one 
full cycle of the protocol lasting 20ms. Fig. [2^ shows the vari- 
ances of the teleported states as a function of gain g for dif- 
ferent input states. The quadratic dependence of the variances 
on g predicted by the model ll23ll fits the experimental data 
very well. It is clear from the figure that for a certain range 
of gain values, the atomic variances are reduced, which indi- 
cates entanglement of the transmitted light with the ensemble 
B |23|. Fig.[2]3 presents the experimental fidelity (blue dots), 
which is above the classical benchmark for n < 7. Instead of 
applying the classical feedback conditioned on the Bell mea- 
surements to the spin B and then performing the verification 
measurement, we can choose to first perform the verification 
measurement on spin B and then apply the classical feedback 
to the results. In theory, those two procedures are equivalent, 
but in the experiment the resulting fidelity for the latter one is 
slightly higher (red dots in Fig.[2j)) since the application of rf 
fields required in the former procedure introduces additional 
technical errors. 

The deterministic teleportation can be used for "strobo- 
scopic" teleportation of spin dynamics from A to B. To il- 
lustrate this specific feature, we performed repeated teleporta- 
tion cycles while varying the amplitude and phase of the input 
state. The results are presented in Fig.|3] The left column dis- 
plays the rf field which is applied to prepare a new spin state 
A each time, after initializing both ensembles to vacuum be- 
tween teleportation runs. The central column shows the read- 
out of the input state evolution of ensemble A and the right 
column shows the read-out of the teleported state evolution. 



The calibration of the input atomic spin state, the Bell mea- 
surement, and the detection of the state teleported onto the 
spin B are performed via polarization homodyning measure- 
ments of the Stokes operator 5*2 = (n+45 — 77,-45) /2 given by 
the difference of photon numbers polarized in ±45° directions 
(Fig[T]) |20||25J. The measured canonical variable for light is 
then defined as 5*2 ~ ■ y where <l> is the driving field 

photon flux, which experimentally means that all measure- 
ments are normalized to shot noise of light. The photocurrent 
is analyzed with a lock-in amplifier at the Larmor frequency 
ft and further computer processed to obtain measurements of 
the temporal modes of interest Vc/s-- AH light pulses for tele- 
portation and read-out always pass through both vapour cells 
(Fig.[T]). For the read out of each individual ensemble the other 
ensemble is detuned from the atom-light interaction by briefly 
detuning the B field in the respective cell. For off -resonant 
light well below saturation used here, the linear transforma- 
tion of light variables after dispersive interaction with atoms 
is given bv l[T9ll25l : 

Vc/s,- =l^-p/x^Cy y'^i^j^ + Cq ' q'^^^ j^ + CAT • F^,^. (1) 

Here, the first term is a contribution of the atomic spin vari- 
able due to Faraday rotation of light polarization, the second 
term is proportional to the input value of the light quadrature 
y of the temporal mode fy and the third term is the contribu- 
tion of the other quadrature of input light q of temporal mode 
fq resulting from back action of light on atoms |23|. All in- 
put light modes are always in a coherent or vacuum state with 
Var{y'^^^ jJ = Var{q^^^ jJ = 1/2. The last term in Eq. [l] 
describes additional noise arising from atomic decoherence 
with Var{Fp^x) = m/2 with m = 1.3 found from the atomic 
spin relaxation |26|. The value of the interaction constant 
is found by calibrating the Faraday rotation caused by the en- 
semble 1 19 |. The constants Cy and Cq are determined by send- 
ing light with displacements of < y^^^^ ^ > and < q^^j^ ^ >, 
storing it in the atomic medium, then reading it out onto an- 
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Figure 3: Teleportation with a dynamically varied input state. 

Left column - rf magnetic field applied to spin A with components 
xa OQ Jy (X V^f^PA (X Jz (X V^f- The center/right columns present 
the read-out of the input/teleported dynamics of the atomic spins A/B 
in vacuum units. Every grey point is a teleportation run with the 
points taken at the rate of ^ 50Hz. The red lines present the running 
average of the grey points. The first/second row is the pa,b/xa,b 
variable and the third row is a two-dimensional plot xa,b,Pa,b- The 
optimal teleportation gain for this evolution is 0.8 which is seen as a 
smaller mean amplitude of the teleported evolution compared to the 
original. 



other pulse < y^^^^ _ > and measuring the ratios. The values 
for Cy, Cq and cat can also be calculated from the model 1231 
based on three experimental parameters: the total trans- 
verse decay rate of the atomic spin state 7, the contribution 
of spin decoherence (spontaneous emission, collisions and in- 
homogeneity of the magnetic field) to this decay time 7extra 
and = (/i + ^)/(m — Z'^ = 6.3 is calculated from 
Clebsch-Gordon coefficients for the atomic transitions and 
experimentally verified |18|. For 5.6mW read out pulses of 
2ms duration and room temperature Cs vapor pressure (ef- 
fective resonant optical depth of 34 for 22mm long cells) 



7 = 99.3 ± 0.2sec~^ and 7extra = 26.3 ± 0.2sec~^ have been 
measured. A unitary contribution to the decay 7 — 7extra is due 
to the collective coupling i^int, which describes the rate of en- 
tanglement generation and the beam splitter interaction 1231 
and depends on the optical depth of the ensemble, the opti- 
cal detuning, and the intensity of the driving field. The mea- 
sured values of n^Cy^Cq agree very well with the predictions 
of the model |23 1 and are n = 0.87, Cy = 0.93, Cq = 0.50 
for our teleportation setting. The last coefficient in the read- 
out equation can be found from the measured parameters as 
cn = Cq - ^J2 • 7extra/(7 " Textra)/^ = 0.17. For the atomic 
state reconstruction the detection efficiencies including opti- 
cal losses riB = 0.80 ± 0.03 and r]A = 0.89 ± 0.03 for ensem- 
bles A/B are taken into account. 
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Supplemental Information 

In the following, we provide additional information on the 
the read-out of the atomic spin state and the protocol used 
for teleportation. In Sec. [A| we discuss the input-output re- 
lations for the interaction of a single cell with light. These 
equations are the basis for the read-out which is used for veri- 
fication. In Sec.[B] we explain the teleportation scheme, state 
the corresponding input-output relations for atoms and light 
and calculate the attainable fidelity. 



relation for these read-out modes are given by 



m 



Cn 



(S.l) 



These equations include atomic decay. The first term on the 
right side is the desired atomic signal. The second and third 
term on the right represent contributions of the light field. 
The first subscript of the photonic operators refers to the fast 
modulation (i.e. modulation with sin(l]t) or cos(ytt)), while 
the second subscript refers to the mode function of the slowly 
varying envelope (as explained below). F^^ and F^^ are 
atomic noise operators, Cy, Cq and cn are real coefficients 
and are given by Eqs. (|S.5|) below. 



The interaction between atoms and light is governed by 
the Hamiltonian H = Ha + i^L + ^int, where Hi^ represents 
the free propagation of the light field along z and Ha is the 
free Hamiltonian of the atomic system. IMore specifically, 
Ha = — ^ {x'^ /2, describes the Larmor precession 

with Larmor frequency Q of the atoms in the magnetic field. 
The interaction Hamiltonian is a combination of a passive 
(beamsplitter-) part H^s, which preserves the number of 
excitations in the system and an active (two mode squeezing) 
part HjMS^ which creates entanglement between atoms and 
light, 



H\r^ 




1 



xy{0) 



(S.2) 



where /i = \{Z ^ u = \{Z — Here, we 

assumed a pointlike atomic ensemble which is located 
at the origin z = 0. The light field is described in 
terms of spatially localized modes y{z^t), q{z^t) with 
[y{z,t),q{z',t)] = c5{z-z') |S1,S2,S3 |. We perform a vari- 
able transformation to a moving frame y{(_^t) = y{ct — (,^t) 
on the spatial argument of the light-field operators (see for 
example IS31 [S4l ). Using this description, integrated light 
modes can be defined by considering an integral over the 
individual pieces of the light pulse with temporally varying 



weighting functions (see Eq. (|S.3|) and Eq. (S.4)). 



A. Light-matter interaction 

We are interested in reading out the spin state of an atomic 
ensemble, which is described in terms of bosonic operators x 
and p. This information is mapped to a coherent light field. 
The interaction of light with atoms, which are rotating in a 
magnetic field with a Larmor frequency ft leads to temporally 
modulated light modes. More specifically, the atomic quadra- 
tures X and p are mapped to sin(l]t) and cos{ftt) modulated 
light modes respectively, which can be accessed individually. 
As explained below, we consider here specific light modes 
with an exponentially falling slowly varying envelope on top 
of the fast sine/cosine modulation. The relevant input-output 



We derive the input-output relations for the light-matter 
interaction given by Eq. ( |S.2| ) in the limit QT ^ 1, where 
T is the total interaction time. Contributions which are on 
the order of {QT)~^ are neglected, which is a very good 
approximation for the experimental parameters considered 
here (compare |[S3l [S5l [56 J ). We include undesired noise 
processes, which lead to a decay of the transverse atomic spin 
at a rate 7extra- In the presence of noise, the optimal slowly 
varying envelope for the read-out is an exponentially falling 
mode e~^\ where 7 = 7s + 7extra- In the experiment, the 
exponentially falling sine and cosine modulated light modes 
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are measured. The corresponding input-output relations are 



given by Eq. ( |S.1| ) above. The Hght and noise modes appear- 
ing in this equation are defined by 



y"j. 




sin(m)\ 



(S.4) 



where fx{t) and fp{t) are atomic noise operators with 
iUt)) = ifpit)) = and {Ut)Ut')) = imUif)) = 
fS{t - t'). Accordingly, var {F'^) = var {F^ = m/2. For 
the interaction Hamiltonian, given by Eq. ( |S.2| ), one obtains 
m = 1.3 lISTllSSll . The modulating functions /y, /q and /n 
read 



fdt) 

with 



2 7s 



1 



V7 VI - e-27T 

"x/TsTextra 1 



V7 VI - e-27T 



-7* 



-7* 



iVN 



Jo 



The coefficients appearing in Eq. ( |S.1| ) are given by 



(S.5) 



It is instructive to consider the limiting case of large » 1, 
with %TZ'^ = const. In this limit, ji = u and we obtain a 
quantum non demolition (QND) interaction with i^int cx pq. 
In the absence of decay, the coefficients of the input-output 



relation Eq. (S.l ) become 

n = Z\f^^ Cy — 1, 



K^/V^, cn = 0. (S.6) 



In this case, the readout equations take a form extensively used 
previously llS3l[S5l 



V3 V </, 



(S.7) 



B. Teleportation scheme 



In this section 
Sec 



we discuss the teleportation scheme. In 
we explain the basic working principle of the pro- 
tocol and provide the corresponding input-output relations. In 



Bl 



Sec B 2 we compute the teleportation fidelity. 



1. Protocol and input-output relations 

A standard teleportation scheme involving the three parties 
Alice, Bob and Charlie consists of the following three steps, 
which allow Alice to teleport a quantum state provided by 
Charlie to Bob. (i) Alice and Bob establish an entangled 
link, which is shared between the two remote parties, (ii) 
Alice performs a Bell measurement on her part of the 
entangled state shared with Bob and an unknown quantum 
state prepared by Charlie, (iii) Alice uses a classical channel 
to communicate the measurement outcome to Bob, who 
performs a local operation on his quantum state conditioned 
on Alice's result. 

In the setup used here (see Fig. la in the main text), the 
quantum state prepared by Charlie is stored in ensemble II. 
This state is teleported to ensemble I, which represents Bob, 
while the light field in x-polarization plays the role of Alice. 
Step (i) in the standard protocol outlined above corresponds 
to the interaction between the light field and the first atomic 
ensemble which results in an entangled state. The distribution 
of entanglement between the two remote sites is realized by 
means of the free propagation of the photonic state. Step 
(ii) corresponds to the interaction of the light field with the 
second ensemble and the subsequent measurement of the 
^/-quadrature of the transmitted light by means of homodyne 
detection. Step (iii) is implemented in the form of a feedback 
operation realizing a conditional displacement on ensemble I 
using radio-frequency magnetic fields. In principle, it is not 
necessary to perform the displacement as long as the outcome 
of the measurement in step (ii) is known. 
The measured values of yl^^_ and yf^l_ are fed back onto the 
first ensemble as explained above. Due to symmetry reasons, 
applying equal gain factors Qx = Qp = g is optimal. This 
yields 



^tele 

,tele ) - Cl ( 1 ) +C„ ( ^jl ) +CN,, ( ^^^^ ] (S.8) 



Pi 



Pi' 




The modes appearing in this equation are defined as in 



Eq. (S.4) with 



/n,I = \/27extra(e' 



-7(T-t) 



-7(2t'-t) 



[l-27s(t'-t)]), 
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/n,ii — 



-^y/TsTextra 9 



-7(T-^) 



-7t _ p-27T^7^^ 



-7* 



(2- 



27s(r-t))]}, 



-7(T-t). 
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The corresponding normalization factors read 



Jo ^ Jo 



and the coefficients appearing in Eq. ( S.8 ( are given by 



ci 



ClI 



v27 7 
), 



27sr , 



27T 



Cn,I = \/^N,I, Cn,II = yA^NJI, Cy = yA/y, Cq = ^ A/q. 

All coefficients and mode functions in the input-output equa- 



tion for the final atomic state Eq. (S.8) carry a bar in order 
to avoid confusion with the coefficients and mode functions 



appearing in the readout equation Eq. ( S.l ) in Sec. A 



2. Teleportation fidelity 

The performance of the protocol is assessed using the aver- 
age fidelity with respect to a Gaussian distribution of coherent 
input states as figure of merit. The fidelity F = \ (^f^^l^j^^) p 
is given by the overlap of the final atomic state in ensemble 
I (Bob's state), , which is described by xf^^ and pf^^ 

and the optimal final state which is defined by the initial state 
in ensemble II (Charlie's state) xjj, pjj. For a given coher- 
ent input state with mean values (xu) and (pu) (and variances 
var(xn) = var(pn) = 1/2 ), the single-shot fidelity is given 
by 



a) 

m 



0.75 

0.7 

0.65 

0.6 

0.55 




Supplementary figure S.l: Average teleportation fidelity F(n) for 
optimal gain g versus width of the distribution of input states n. The 
lowest (red) line in both panels represents the classical benchmark, 
a) Teleportation fidelity for the measured experimental parameters 
(see methods), b) Maximum teleportation fidelity in the absence of 
losses. The curve in the middle (black) shows the attainable telepor- 
tation fidelity for the interaction used in the experiment (Z = 2.5). 
The upmost curve (blue) depicts the QND-fidelity (for Z ^ oo, 
j^TZ^ — const) for exponentially shaped driving pulses (compare 
text, Sec. [62]). 



such that the average fidelity F{n) with respect to a Gaussian 
distribution with width n. 



F{n) 



1 r^r 



d{xii)d (pii )F{{xii), {pii )) e ' 
V2 



. lfqi)_t(pn)_ 



l + 2var (xfi^)+2n(^l 
V2 



(S.9) 



{xii) 



l + 2var(pf«)+2n(l-|gtLp 



Both ensembles are initialized in a coherent spin state with 
{xi) = (xii) = and (x^) = (x^) = 1/2. The photonic 
modes are also initially in the vacuum state, such that 



\{xi)/{xii)\ = \{pi)/{pii) \ = Cii, 
var(xfi^) = =^(ci' ' 
var(pfi^) 



^N,II 



F{{xu)Apii)) 



var(xr) 

Fig. S.l shows the average teleportation fidelity F{n) in 
comparison to the classical limit Fdas = (1 + n) / (1 + 2n), 
which cannot be surpassed by classical means fS9l lSlQ|[Slll . 
This figure also displays the average fidelity which can be 
achieved using a QND-interaction (i.e. for very large detun- 
(l + 2var(xf^^)) (l + 2var(pf^^)) i^g, see SecfAl if the classical driving pulses are modulated 



-(l(^ii)l-K^^'^)l)^ 

g l + 2var(xfle) 



-(KPll)l-l(pf'^)l)^ 
l + 2var(ptele) 
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in time. These results have been obtained by considering dif- teraction. The fideUty is optimized with respect to / and /„. 
ferent exponential functions fi{t) oc e^^^ and /ii(t) oc e-^"^ for Fig. S.lb shows that fidelities close to one can be obtained in 
the pulse shape of the classical field in the first and second in- principle. 
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